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Abstract. We study dynamic behavior of Potts model with invisible states near the first-order 
phase transition temperature. We focus on melting process starting from the perfect ordered 
state. This model is regarded as a standard model to analyse nature of phase transition. 
We can control the energy barrier between the ordered state and paramagnetic state without 
changing the symmetry which breaks at the transition point. We calculate time-dependency 
of the order parameter, density of invisible state, and internal energy. They show two-step 
relaxation behavior. We also consider the relation between the characteristic melting time and 
characteristic scale of the energy barrier by changing the number of invisible states. We find 
that characteristic melting time increases as the energy barrier enlarges in this model. Thus, 
this model is regarded as a fundamental model to analyze dynamic behavior near the first-order 
phase transition point. 



1. Introduction 

Frustration causes many interesting static and dynamic behavior which are not observed in 
unfrustrated systems because of pecuhar density of states [ItilU]. In two-dimensional frustrated 
systems, there have been found many nontrivial phase transitions such as order by disorder 
[TT1[T2], reentrant phase transition [T^TT], topological phase transition [TH], and novel type of 
first-order phase transition. Recently, strange first-order phase transitions have been found in 
two-dimensional frustrated continuous spin systems |19H22j . 

In [TU], the authors studied equilibrium properties of the classical Heisenberg model on 
triangular lattice with nearest neighbor ferromagnetic interaction Ji and third-nearest neighbor 
antiferromagnetic interaction J3. They found that a first-order phase transition with threefold 
symmetry breaking occurs at finite temperature. This looks a strange phase transition, since 
phase transition with threefold symmetry breaking is often second-order phase transition on 
two-dimensional lattice e.g. the three-state ferromagnetic Potts model or the three-state 
ferromagnetic clock model [23j. After this study, similar nature of first-order phase transition 
have been found by a number of researchers [20H22j . Stoudenmire et al. found a first-order phase 
transition with breaking of threefold symmetry in Ji — J3 model with biquadratic interaction on 
triangular lattice. Okumura et al. also found similar first-order phase transition in frustrated 
Ji — J2 model on hexagonal lattice. Such a first-order phase transition sometimes takes place 
in two-dimensional frustrated systems with a number of competed interactions. It is an open 



problem why such a first-order phase transision appears in some two-dimensional frustrated 
systems. 

To consider this problem, we constructed a model which exhibits a first-order phase transition 
with threefold symmetry breaking by introducing new kind of parameter into the standard 
ferromagnetic three-state Potts model [MIES]. We introduced the invisible states which does 
not contribute to the internal energy. The ground state of the Potts model with invisible state 
and its degeneracy are the same as that of the standard Potts model. We found that the first- 
order phase transition with threefold symmetry breaking occurs by just adding the invisible 
states into the standard three-state ferromagnetic Potts model. The standard Potts model 
is regarded as a fundamental model to consider the properties of phase transition in statistical 
physics [23p26| . It is believed that the Potts model with invisible states is a potential model which 
can clarify inherent nature of first-order phase transition with breaking of threefold symmetry 
on two-dimensional lattice. 

Our aim of the present study is to clarify dynamic properties of Potts model with invisible 
states. It is interesting topic in statistical physics to study dynamical nature of the systems 
which exhibit a first-order phase transition, since there are some similar points between these 
systems and glassy systems [271[28]. It is expected that some dynamic nature in glassy systems 
can be explained in terms of dynamic behavior of the systems which have a first-order phase 
transition. Potts model with invisible states can be considered a potential model to study 
dynamical nature of first-order phase transition. This is because we can control the energy 
barrier between the ordered state and paramagnetic state without changing the symmetry which 
breaks at the transition temperature by modulating the number of invisible states. Especially, 
we focus on melting process in this model. In this paper, we study dynamics of the Potts model 
with invisible states at fixed temperatures which are above the transition temperature. 

2. Model 

We consider the ferromagnetic Potts model with invisible states on square lattice. The 
Hamiltonian of this model is 

(jj) Q=l 

at = 1,- ■ ■ ,q,q + l,- ■ ■ ,q + r, (2) 

where denotes the nearest neighbor pairs on square lattice. In this paper we take J as an 
energy unit. If and only if 1 < Cj = aj < q, the interaction works. Obviously this model for r = 
corresponds to the standard g-state ferromagnetic Potts model. We call the state in 1 < ai < q 
"colored state" and the state m q + 1 < < q + r "invisible state" from now on. Hereafter we 
call this model (g,r)-state Potts model. The invisible state does not affect the internal energy. 
It should be noted that the number of ground states of the (g,r)-state Potts model is the same 
as that of the standard ferromagnetic g-state Potts model. Then the (g,r)-state Potts model 
exhibits a phase transition with q-fold symmetry breaking. To change the number of invisible 
states r corresponds to changing form of the density of states. The invisible states contributes 
the entropy as logr. Then, it is expected that the order of phase transition can be changed 
by adding the invisible states. Actually, a first-order phase transition with q-fold symmetry 
breaking occurs for large enough the number of invisible states r even for q = 2,3, and 4 on 
two-dimensional lattice [23]. As the number of the invisible states r increases, the latent heat 
increases and the transition temperature decreases [25] . 



3. Result 

In this paper, we study dynamics of the {q,r)-state Potts model on square lattice whose size is 
A'^ = 128 X 128 at fixed temperatures which are above the transition temperature following [28j. 
We impose periodic boundary condition. We prepare independent 1024 samples for obtaining 
data with high accuracy. The initial state is set to be perfect ordered state such as ctj = 1 for all 
i. We adopt single-spin- flip Metropolis type of Monte Carlo method as the time-evolution rule. 

To consider dynamic behavior of the (g,r)-state Potts model, we define the time-dependent 
order parameter m{t), density of invisible states pinvit), and internal energy e{t) as follows: 
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where t denotes Monte Carlo step. Figure [1] shows the time-evolution of the order parameter 
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Figure 1. Time-evolution of the order parameter density of the invisible states (pinv(i))) 

and the internal energy {e{t)) of the (3,27)-state Potts model starting from the completely 
ordered state ai{0) = 1 (for all i) for fixed different temperatures. From left to right, the 
temperatures are T = + 0.5, + 0.2, + 0.1, + 0.05, + 0.02, + 0.01, + 0.005, 
Tc + 0.002, Tc + 0.001, Tc + 0.0005, + 0.0002, and + 0.0001, where = 0.58513 t25j. 
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Figure 2. The dynamical susceptibility of the order parameter Xmit) and density of invisible 
states Xpinv(0 of (3,27)-state Potts model starting from the completely ordered state 
(Ji(0) = 1 (for all i) for fixed different temperatures. From left to right, the temperatures 
are T = Tc + 0.5, Tc + 0.2, + 0.1, + 0.05, + 0.02, + 0.01, Tc + 0.005, + 0.002, 
Tc + 0.001, Tc + 0.0005, Tc + 0.0002, and + 0.0001, where = 0.58513 [25]. 



{m{t)), density of the invisible states {pinvit)), and the internal energy (e(t)) of the (3,27)- 
state Potts model starting from the completely ordered state (Tj(0) = 1 (for all i) for fixed 
different temperatures. Here (•) denotes ensemble average. In Fig.[Tl there are obvious two-step 
relaxations. Two-step relaxation often appears in systems where a first-order phase transition 
takes place and also in glassy systems |29] . Melting occurs after stabilization in the ordered state 
which corresponds to plateau region. As the temperature approaches to the transition point, 
the plateau region enlarges and as a result, the characteristic melting time becomes long. 

Next we consider the dynamical susceptibility to consider the characteristic melting time 
systematically. The dynamical susceptibility of the time-dependent physical quantity A{t) is 
defined as 

XA{t) = Np{{A{tf)-{A{t)f), (6) 

where /3 denotes the inverse temperature. We calculate the dynamical susceptibility of the order 
parameter Xm (t) and density of the invisible states Xpinv (^) • Figure [2] shows the dynamical 
susceptibility Xm(t) and Xpinv(^) starting from the completely ordered cTj(O) = 1 (for all i) 
state for fixed different temperatures of the (3,27)-state Potts model. In Fig. [21 the dynamical 
susceptibilities have a peak which indicates characteristic melting time. Characteristic melting 
time defined from Xm{t) and Xpinv(*) almost same. Then we define Tmax as the peak position 
of dynamical susceptibility of the order parameter Xmit)- As the temperature decreases, the 
peak height increases. This result means that characteristic length scale becomes large as the 
temperature approaches to the transition point, since the peak height relates to characteristic 
length scale. The characteristic time Tmax also increases as the temperature decreases. These 
dynamical susceptibilities are quantitatively similar with ^4 which is often used in analysis of 
glassy systems [28] . 

We also study temperature-dependency of Tmax for the (3,25)-state Potts model and (3,27)- 
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Figure 3. The peak position of the dynamical susceptibihty of the order parameter Tmax as a 
function of 1/(T — T^), where Tc is the transition temperature. The red squares and blue circles 
indicate the case for the (3,25)-state Potts model and the (3,27)-state Potts model, respectively. 
Transition temperatures for the (3,25)-state Potts model and the (3,27)-state Potts model are 
0.59630 and 0.58513, respectively [25]. The dotted line indicates 1/(T - Tc)^ ^. 
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state Potts model to consider the effect of the number of invisible states. Figure [3] shows the peak 
position as a function of 1/(T — Tc), where Tq denotes transition temperature. The transition 
temperatures for (g,r)=(3,25) and (3,27) are 0.59630 and 0.58513, respectively [25]. In high 
temperature region, the peak position obeys Tmax oc 1/(T — Tc)^'^ whereas near the transition 
temperature, the characteristic time Tmax becomes longer rather than 1/(T — Tc)^'^. In high 
temperature region, behavior of Tmax is the same both for (3,25)-state Potts model and for 
(3,27)-state Potts model. However, there is difference near the transition temperature. The 
characteristic time Tmax for the (3,27)-state Potts model is larger than that for the (3,25)-state 
Potts model, since the energy barrier between the ordered state and paramagnetic state for the 
(3,27)-state Potts model is larger than that for the (3,25)-state Potts model |24ll25j . From this 
result, we expect that Tmax enlarges as the number of invisible state r increases, since the latent 
heat increases as r increases [25j. Thus, the Potts model with invisible states can be regarded 
as a standard model to analyze dynamic behavior of systems which exhibits a first-order phase 
transition. 

4. Conclusion and future perspective 

We studied dynamic properties of the Potts model with invisible states at fixed temperatures 
which are above the transition temperature. In this paper, we focused on melting process which 
is an important nature of first-order phase transitions. First, we considered the dynamical 
properties of the order parameter, density of invisible states, and the internal energy. We found 
that there is two-step relaxation which is typical behavior in systems which exhibit a first-order 
phase transition and also glassy systems. Next we calculated the dynamical susceptibility of the 
order parameter and density of invisible states. As the temperature approaches the transition 
temperature, the peak height of the dynamical susceptibility grows. Furthermore, we studied 
the relation between the peak position of the dynamical susceptibility of the order parameter 
Tjnax and temperature. We found that the peak position obeys Tmax oc 1/{T — Tc)^'^ in high 
temperature region, whereas this relation breaks near the transition temperature. The growth 
speed of Tmax depends on the energy barrier between the ordered state and paramagnetic state. 
In the Potts model with invisible states, we can control the energy barrier between the ordered 



state and paramagnetic state without changing the symmetry which breaks at the transition 
point. We should study the relation between the energy barrier and the time-scale and length- 
scale in this model more carefully. Mori and Nakada have studied dynamical properties in this 
model systematically to clarify whether there is new universal behavior near the first-order phase 
transition point |30j . 

It is also important topic to clarify whether the Potts model with invisible states relates to 
strange type of first-order phase transition in two-dimensional frustrated systems. It will be 
reported elsewhere. 
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